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I .  INTRODUCTION 


This  working  paper  addresses  the  possibility  of  directly 
inverting  teleseismic  ground  motion  in  order  to  obtain  explo- 
sion yield.  We  propose  to  accomplish  this  inversion  via  a 
deterministic  model  which  predicts  the  teleseismic  ground 
motion  for  a nuclear  explosion. 

A schematic  of  the  model  is  shown  in  Figure  1.1.  It 
is  the  explosion  yield,  which  controls  the  amplitude  of  the 
equivalent  elastic  source,  that  we  wish  to  obtain  based  on 
a comparison  between  the  predicted  and  observed  ground  motion. 

Obviously  this  approach  is  nonsense  unless  the  remain- 
ing key  parameters  in  the  model  are  known.  For  a given  re- 
ceiver structure  and  instrument  response  the  key  parameters 
are  as  follows: 

1.  The  material  properties  and  constitutive  model 
of  the  near  source  rock  environment. 

2.  The  pP  delay  times. 

3.  The  anelastic  properties  of  the  earth's  upper 
mantle . 

What  can  we  learn  from  the  data  exchange?  Our  intent 
is  to  use  the  data  to  choose  an  analog  rock  whose  material 
properties  are  consistent  with  the  data  furnished  in  the  ex- 
change. This  provides  enough  information  to  permit  an  a 
priori  calculation  of  the  equivalent  source. 

If  we  can  determine  the  equivalent  source  for  the  near 
source  rock  environment,  then  the  pP  delay  time  and  the  an- 
elastic properties  of  the  earth's  upper  mantle  required  Y ' 
our  model  may  be  found  from  the  details  of  the  recorded  ground 
motion. 
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Figure  1.1.  Schematic  representation  for  the  deterministic  prediction  of  the 
teleseismic  ground  motion  from  a nuclear  explosion. 


II.  THE  EQUIVALENT  ELASTIC  SOURCE 


Figures  2.1  and  2.2  show  the  calculated  far-field  dis- 
placement spectra  of  the  indicated  rock  types.  In  all  cases 
the  device  yield  was  150  kt  and  the  DOB  equaled  1860  ft  (350 
. Table  1 lists  the  material  properties  used  for  each 
calculation.  The  shot  name  for  the  event  located  in  each 
material  is  also  given  in  the  table.  We  should  note  that 
measured  material  properties  were  used  in  the  source  calcula- 
tions. These  sources, when  used  in  conjunction  with  the 
model  in  Figure  1.1,  explain  both  the  free-field  and  tele- 
seismic  data  available  to  us. 

Since  there  is  a factor  of  ten  difference  in  amplitude 
at  1 Hz  between  the  maximum  and  minimum  source  functions,  it 
is  important  to  subdivide  these  sources  into  classes.  Table 
2 performs  this  separation.  Figure  2.3  shows  the  source 
functions  for  all  rock  types  except  Class  I (207)  and  illus- 
trates the  natural  decomposition  of  sources  into  their  respec- 
tive classes.  It  is  our  intention  to  use  the  information 
from  the  data  exchange  to: 

1.  Classify  the  near  source  rock  environment. 

2.  Perturb  the  material  properties  in  the  class 
to  fit  the  given  depth  of  burial,  density, 

P wave  velocity  and,  hopefully,  S wave  velocity. 

3.  Calculate  the  range  of  equivalent  sources  which 
are  possible  in  this  class.  This  range  will 
include  reasonable  variations  in  material  strength 
and  air-filled  porosity  consistent  with  past 
experience  from  similar  rock  types. 
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Figure  2.3.  Source  functions  for  all  rock  types  except  207 


TABLE  2 

I 


Calculation 

Class 

Description 

Rock  Type 

a 

V 

¥ (1.0) 

T 

207 

I 

Weak,  saturated 

Tuff 

2.4 

40 

19 

201 

II 

Strong,  saturated 

Granite 

5.3 

207 

12.0 

203 

II 

Strong,  saturated 

Rhyolite 

4.2 

170 

15.3 

204 

11 

Strong,  saturated 

Tuff 

3.1 

90 

11.2 

206 

III 

Weak,  dry 

Shale 

3.9 

150 

4.6 

202 

IV 

Strong,  dry 

Granite 

5.3 

207 

2.5 

205 

IV 

Strong,  dry 

Sandstone 

4.3 

204 

2.0 
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III.  TELESEISMIC  AMPLITUDE  DEPENDENCE  ON  THE  SOURCE 


The  fundamental  output  of  our  nonlinear  finite  dif- 
ference calculations  is  the  reduced  displacement  potential 
(4'(t))  representation  of  the  explosion  source.  In  Figure 

A 

2.1  we  plotted  |'t'(oj)|  for  seven  150  kt  explosion  detonated 
in  a representative  range  of  emplacement  materials.  These 
source  spectra  span  a considerable  range  of  amplitudes  and 
spectral  shapes.  How  do  the  source  spectra  translate  into 
short  period  seismogram  amplitudes  at  teleseismic  ranges; 
that  is,  into  m^?  The  relation  between  the  two  is  discussed 
in  this  section. 

In  a homogeneous  space,  the  displacement  is  related 
to  the  reduced  displacement  potential  in  the  frequency  do- 
main by 


A 

= (! 


(<j)  , Y (o>) 
Ra  „2 


-ikaR 
e , 


(3.1) 


where  k = to/a.  In  the  far-field  the  R~ 1 term  dominates, 
a 

The  vertical  component  of  the  teleseismic  ground  motion  is 
then  related  to  the  source  region  displacement  through  a 
series  of  transfer  functions.  The  relation  may  be  written 


uR  = us  * Ts(u)  * TR(u)  • Tm(6j)  • I(w),  (3.2) 


where 

Ts(w) 

TrU) 

V“> 

I(u>) 


the  influence  of  the  crust  in  the  source  region, 
the  influence  of  the  crust  in  the  receiver  region, 
the  influence  of  the  upper  mantle, 
the  seismograph  response. 
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The  transfer  functions  representing  the  effect  of  the 
earth's  crust  and  upper  mantle  may  be  quite  complex.  Here 
we  shall  isolate  the  source  effect  by  considering  only  the 
most  elementary  forms  for  these  transfer  functions. 

Assume  that  the  upper  mantle  transfer  function  takes 
the  following  form: 

V“>  = tgs  e"  f"T/Q  <3-3> 


where  TGg  is  the  effective  1/R  or  geometric  spreading  fac- 
tor for  the  upper  mantle  path.  The  T/Q  is  taken  to  be  a 
constant  characterising  the  path.*  For  the  TR (co ) we  simply 
use  a constant  (a  ) which  converts  the  emerging  body  waves 
to  vertical  motion  at  the  free  surface.  Then  (3.2)  reduces 
to 


A 


aR  TGS  I(U))  TSU) 


$ (w)  -irf  T/Q 

-z — e ' 


(3.4) 


and  synthetic  seismograms  may  be  computed  by  transforming 
uR  to  the  time  domain. 

We  now  specify  the  common  factors  (a_,  T__,  I (oj)  , T/Q) 
and  compute  synthetic  seismograms  for  the  range  of  source  func- 
tions represented  in  Figure  2.1.  The  relative  amplitude  and 
frequency  content  of  these  seismograms  will  be  dependent  only 
on  the  character  of  the  source  and  source  crustal  structure. 
Seismograms  are  computed  for  four  cases  of  increasing 
complexity : 


A complex  exponential  representing  the  dispersion  due  to  Q 
(Strick,  1970)  is  also  included  in  the  calculations  but  is 
left  out  of  the  formula  for  clarity. 
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1. 


The  source  region  crust  is  a homogeneous  material 
of  the  source  material.  That  is,  T c (u>)  = 1. 

u 

2.  A different  layered  earth  model  represents  the 
crust  for  each  source.  The  top  layer  is  a 

half space  of  the  source  material.  That  is,  there 
is  no  free  surface . A few  kilometers  below  the 
source  all  crustal  models  are  the  same.  Between 
the  source  and  this  depth  a few  layers  are  in- 
serted to  represent  a fairly  smooth  increase  in 
velocity  with  depth. 

3.  The  same  structures  are  used  as  in  2 but  now  a 
free  surface  is  introduced  in  such  a way  as  to 
fix  the  delay  time  between  the  direct  (P)  and 
free  surface  reflected  (pP)  waves.  That  is, 

the  burial  depth  is  inversely  proportional  to  the 
a for  the  source  material.  (A  small  correction 
is  required  to  account  for  differences  in  the 
takeoff  angle.) 

4.  Calculations  like  3 are  carried  out  but  now  with 
depth  held  constant.  Then  the  P-pP  lag  time  is 
inversely  proportional  to  a. 

The  method  used  to  compute  Tg  (w)  for  the  latter  two 
cases  is  essentially  identical  to  that  developed  by  Fuchs 
(1966) . This  method  has  been  generalized  to  more  complex 
sources  by  Bache  and  Harkrider  (1976)  and  the  actual  imple- 
mentation is  as  outlined  there.  For  the  first  two  cases 
listed  above,  the  Haskell  matrix  formulation  was  modified 
to  allow  the  first  layer  to  be  a half space  rather  than  having 
a free  surface  (unpublished  work) . 

The  common  factors  for  the  synthetic  seismogram  cal- 
culations are  as  follows: 

aR  = 1*7' 
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T = 0.95  x 10“ 4 (a  typical  value  for  distances  like 

4000  km) , 

T/Q  = 0.7  (perhaps  a typical  value  for  Eurasian 
travel  paths) . 

For  the  instrument  we  took  response  curves  provided  by  the 
VSC  Project  Officer,  but  these  are  not  substantially  dif- 
ferent from  LRSM  nominal  curves.  All  calculations  are  done 
at  a yield  of  150  kt. 

The  seismograms  for  cases  1 and  2 are  shown  in  Figure 
3.1.  Case  1 is  the  simplest.  Essentially,  we  are  taking 
the  far-field  term  from  the  whole  space  equation  (3.1)  and 
transforming  to  the  time  domain.  We  also  include  the  modifi- 
cation by  a realistic  Q and  a typical  instrument.  Thus,  the 
seismograms  in  Figure  3.1a  represent  the  direct  P wave, 
unmodified  by  any  layering  effect.  Parenthetically,  we 
point  out  that  the  source  displacement  was  computed  at 
= 20  km  from  the  source,  accounting  for  the  travel  time  dif- 
ferences apparent  in  the  figure. 

The  source  material  is  different  for  each  of  the 
seismograms  of  Figure  3.1a.  While  the  material  in  the  source 
vicinity  may  vary  over  the  range  of  materials  represented, 
we  suppose  that  at  some  depth  the  underlying  rocks  are  pretty 
much  the  same.  The  direct  P wave  will  certainly  be  affected 
by  variations  in  the  impedence  mismatch  between  the  source 
material  and  the  underlying  common  geology. 

In  computing  the  seismograms  of  Figure  3.1b  we  assumed 
a different  crustal  model  for  each  source.  These  models  are 
tabulated  in  Table  3.  We  see  that  these  are  very  simple 
crustal  models  with  the  source  material  occurring  in  a top 
layer  1.3  km  thick.  This  is  underlain  by  crustal  granite. 

For  all  the  calculations  a horizontal  phase  velocity, 
c = 12.69  km/sec,  was  assumed.  This  translates  into  a take- 
off angle  of  28.2°  at  the  base  of  the  crust.  The  takeoff 
angle  at  the  source  is  found  by  applying  Snells'  law;  that  is, 
0 = sin”1  a/c. 
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yield  is  150  kt.  There  is  no  free  surface  in  the  calculation.  The  per 
of  the  first  (b)  and  second  (d)  cycles  on  the  record  is  indicated.  The 
numbers  at  the  left  on  these  and  subsequent  synthetic  seismograms  are 
amplitudes  in  microns  at  one  Hertz. 


♦ 


In  summary,  the  seismograms  of  Figure  3.1  are  for  the 
direct  P wave;  first/,  including  no  layering  effects  and  second, 
including  the  propagation  from  the  source  material  into  the 
underlying  common  geology.  From  each  of  the  seismograms  we 
make  amplitude  measurements  in  the  standard  way;  that  is, 
we  measure  a peak-tj^-peak  amplitude  and  time  for  both  the 
b and  d phases.  The  periods  T^  and  T^  are  then  twice  the 
peak-to-peak  time  and  the  amplitude  is  corrected  for  the 
instrument  response  at  this  period. 


How  do  the  amplitudes  scale  with  the  source  material? 
For  case  1 it  is  immediately  clear  that  each  frequency  com- 

/\  A 

ponent  of  uR(co)  scales  with  | $ (u)  | /a . All  other  factors 
are  held  constant  in  the  calculation.  The  only  question 
is  then  how  much  the  phase  of  the  spectrum,  the  Q effect, 
the  instrument  and  the  measuring  procedure  itself  obscure 
this  scaling  relationship. 


When  layering  is  introduced,  the  scaling  is  not  so 
obvious.  Bache,  et  al . (1975,  Appendix  A)  show  that  in  this 

A . "• 

case  uR(co)  scales  with  a | V (a>)  | . In  brief,  the  argument  is 
that  the  far-field  term  in  the  whole  space  equation  (3.1) 
is  replaced  by 


A [I  (a;)  1 
L a 


(3.5) 


where  ug  (a>)  is  now  the  displacement  spectrum  in  the  common 
material  underlying  all  the  source  materials.  The  A is  the 
reflection  coefficients  and  L replaces  R and  represents  the 
spreading.  Bache,  et  al . (1975)  show  that  A = ot/ab  and 
L ~ ab/a  where  is  the  P wave  velocity  of  the  common 
material.  It  then  follows  that  u_  (<d)  , and  therefore  u_  (w)  , 

£ b K 

scale  with  a | (a> ) | . 

From  the  T^  and  T^  values  given  in  Figure  3.1  we  see 
that  the  dominant  frequency  for  both  the  b and  d phases  is 
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TABLE  3 

CHJSTAL  MODELS  FOR  SOURCES  201-206 


Layer 

Depth  (km)  Thickness  (km) 

a (km/sec) 

8 (km/sec) 

p (q/cm  ) 

: kkk 

Sources  201  - 202 

*** 

l 

1.3 

1.3 

5.33 

2.78 

2.67 

2 

4.1 

2.8 

5.80 

3.45 

2.80 

3 

20.0 

15.9 

6.0 

3.50 

2.80 

*** 

Source  203  *** 

1 

1.3 

1.3 

4.20 

2.63 

2.45 

2 

2.0 

0.7 

4.7 

2.70 

2.60 

3 

2.7 

0.7 

5.4 

2.80 

2.70 

4 

4.0 

1.3 

5.8 

3.45 

2.80 

5 

20.0 

16.0 

6.0 

3.50 

2.80 

•kick 

Source  204  *** 

i 

1.3 

1.3 

3.11 

2.02 

2.2 

2-5 

As  in  crust  for 

Source  203 

*** 

Source  205  *** 

1 

1.3 

1.3 

4.28 

2.85 

2.51 

2-5 

As  in  crust  for 

Source  203 

*** 

Source  206  *** 

1 

1.3 

1.3 

3.87 

2.37 

2.61 

2-5 

As  in  crust  far 

Source  203 

15 

L 


near  1.4  Hz.  The  accuracy  of  the  scaling  relationship  dis- 
cussed above  can  then  be  assessed  by  plotting  the  b and  d 

i * i 

amplitudes  versus  | ¥ (f=l . 4) | /a  for  case  1,  where  the  source 
region  is  homogeneous  and  versus  aj$(f-1.4)|  for  case  2 
where  a layered  geology  is  introduced  below  the  source. 

The  plots  are  shown  in  Figure  3.2. 

The  average  ratio  of  abscissa  to  ordinate  was  computed 
for  each  of  the  data  sets  of  Figure  3.2  and  a line  represent- 
ing the  average  slope  appears  on  the  plot.  An  estimate  of 
the  deviation  from  this  line  was  made  and  the  maximum  devia- 
tion (e  ) is  also  indicated  on  the  plot. 

max 

Figure  3.2a  shows  that  the  scaling  of  source  displace- 

A 

ment  spectral  amplitudes  with  provides  an  excellent 

approximation  to  the  scaling  of  the  time  domain  amplitudes 
from  the  synthetic  seismograms  of  Figure  2.1.  When  layering 
is  introduced,  we  see  in  Figure  3.2b  that  the  spectral 

A 

scaling  with  a|*F(u))|  provides  just  as  good  an  approximation. 

In  summary.  Figure  3.2b  provides  assurance  that  we  under- 
stand how  the  direct  P wave  scales  with  the  source  spectrum 
and  elastic  properties  of  the  source  material.  Further, 

Figure  3.2a  shows  the  kinds  of  errors  that  are  inherent  in 
the  procedure,  even  when  we  know  exactly  how  each  frequency 
component  of  the  seismogram  scales. 

The  next  complication  to  be  introduced  is  the  free 
surface  effect.  Inclusion  of  the  free  surface  in  our  calcula- 
tions actually  has  two  effects.  Most  important  is  the  addi- 
tion of  the  large  negative  phase,  pP,  to  the  seismogram. 

Also,  the  energy  reflected  by  the  impedence  contrast  between 
the  source  and  basement  materials  then  reflects  back  from  the 
free  surface  and  can  influence  the  seismogram's  character. 

Tbe  importance  of  this  latter  effect  depends,  of  course,  on 
the  severity  of  the  impedence  mismatch  and  on  its  depth.  In 
general,  though,  the  effect  is  rather  small. 


16 


f 


The  seismograms  for  cases  3 and  4 , where  the  free  sur- 
face is  included  in  the  calculation,  are  shown  in  Figure  3.3. 
For  computing  these  seismograms  the  crustal  structures  of 
Table  3 were  used..  First,  in  Figure  3.3a  we  take  the  delay 
time  between  P and  pP  (the  P-pP  lag)  to  be  constant  at  0.6 
seconds.  Second,  in  Figure  3.3b  we  consider  the  depth  to  be 
fixed,  thus  the  P-pP  lag  varies  with  a.  Since  the  source 
material  is  assumed  to  extend  to  the  free  surface  (Table 
3)  , 

PH 

P - pP  lag  = cos0  , (3.6) 

where  H is  source  depth  and  0 = sin-1  a/c,  where  c = 12.69 
km/sec  for  these  calculations.  We  should  point  out  that  to 
obtain  the  desired  values  of  0.6  seconds  for  the  lag,  some 
of  the  crustal  structures  of  Table  3 were  modified  by  in- 
creasing the  thickness  of  the  first  layer  while  leaving  the 
rest  of  the  structure  unchanged. 

Measurements  of  the  amplitude  and  period  of  the  b and 
d phase  were  made  for  the  seismograms  of  Figure  3.3,  with 
the  exception  that  no  consistent  d phase  could  be  measured 
on  the  Source  204  record  of  Figure  3.3b.  We  saw  in  Figure 
3.2  that  the  scaling  of  the  direct  P wave  is  quite  accurately 

A 

represented  by  ot  | 'F  ( f ) | where  f is  the  dominant  frequency 

P P 

in  the  transient  signal.  Let  us  then  see  how  much  this 
scaling  is  compromised  by  the  addition  of  the  free  surface 
effect,  especially  pP.  We  note  that  the  composite  pulse 
making  up  the  d phase  has  a period  quite  different  than  the 
0 . 7 second  periods  that  we  know  dominate  the  P and  pP 
phases  taken  separately. 

A 

The  b and  d amplitudes  are  plotted  versus  a | ¥ (f =1 . 4 ) | 

in  Figure  3.4.  The  plot  format,  procedure  for  fitting  the 

line  and  determining  £ are  as  in  Figure  3.2.  First,  con- 

max 

sider  Figure  3.4a,  the  case  where  the  p-pP  lag  is  constant. 
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We  note  that  the  agreement  of  the  b phase  data  with  the  scal- 
ing relation  is  about  as  good  as  in  Figure  3.2.  For  the 

0.6  second  lag  time  the  pP  phase  has  a small  and  consistent 
effect  on  the  b phase.  The  agreement  begins  to  deteriorate 
for  the  d amplitude,  but  is  still  fairly  good.  This  is  an 
indication  that  the  pP  effect  is  consistent  from  event  to 
event.  We  should  mention  that  scaling  with  the  source  spec- 
trum at  a frequency  closer  to  that  of  the  d period,  say 
f = 1.0  Hz,  gives  worse  agreement. 

For  the  constant  depth,  variable  lag  time  case,  the 

A 

data  is  plotted  against  a | $ (f=l . 4 ) | in  Figure  3.4b.  Once 
again,  for  the  b phase  the  agreement  is  reasonably  good 
because  the  pP  arrives  too  late  to  have  a substantial  ef- 
fect. For  the  d phase,  however,  agreement  with  the  scaling 
relation  can  be  qu_te  poor.  For  example,  compare  sources 
201  and  203  or  202  and  205.  Clearly,  the  pP  phase  has  a 
major  effect  on  the  d amplitude. 

In  this  section  all  features  of  the  source-receiver 
travel  path  have  been  held  fixed  except  for  the  crustal 
structure  in  the  vicinity  of  the  source.  We  saw  that  the 
scaling  of  seismogram  amplitudes  is  given  to  quite  good 

A 

accuracy  by  ot  | 'F  ( f ) | except  when  this  scaling  is  severely 
compromised  by  variations  in  the  P-pP  lag  time.  In  principle 
we  should  be  able  to  correct  for  this  latter  effect  and  this 
is  the  subject  of  the  next  section. 

SUMMARY 

1.  Given  the  equivalent  source,  then  teleseismic 
ground  motion  is  directly  proportional  to  a. 

The  single  most  important  near  source  material 
property  is  a.  We  should  insist  that  this 
parameter  be  measured  in  situ. 

2.  If  possible,  use  the  b amplitude  to  obtain  magni- 
tude. For  the  yields  and  DOB  of  interest,  this 
amplitude  is  not  significantly  contaminated  by 
pP. 
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IV.  BODY  WAVE  AMPLITUDE  CORRECTION  FOR  pP 

In  the  previous  section  we  saw  that  the  free  surface 
reflected  phase,  pP,  has  a significant  effect  on  measured 
body  wave  amplitudes.  This  effect  is  investigated  systemati- 
cally in  this  section.  We  compute  a series  of  theoretical 
seismograms  for  which  all  parameters  are  fixed  except  the  depth 
of  burial  or  P-pP  lag.  The  source  is  chosen  to  be  that  which 
was  computed  specifically  for  the  Pahute  Mesa  rhyolite  event 
MAST  (Figures  2.1  and  2.2).  The  source  region  c.'ustal  model 
is  given  in  Table  4. 

The  synthetic  seismogram  computations  were  carried 
out  in  a similar  fashion  to  those  of  the  previous  section. 
However,  a more  complex  procedure  was  used  for  specifying 
the  upper  mantle  and  receiver  crust  transfer  functions, 

Tm(oj)  and  Tr(oj)  • A more  detailed  description  of  this  pro- 
cedure has  been  given  in  previous  S^  reports  (e.g.,  Bache, 
et  al . , 1975;  Barker,  et  al . , 1976;  Bache,  et  al ■ , 1976). 
Haskell  matrix  methods  (Haskell,  1962)  are  used  for  comput- 
ing Tr(w).  However,  the  receiver  structure  was  taken  to  be 
an  average  continental  crust  which  has  almost  no  effect  on 
frequency  content  of  the  seismogram.  For  the  upper  mantle 
elastic  response  a generalized  ray  theory  program  is  used 
(Wiggins  and  Helmberger,  1974).  The  upper  mantle  model  was 
a slightly  modified  version  of  C2  given  by  Anderson  and 
Hart  (1976)  and  the  epicentral  distance  was  3730  km.  At 
this  range  the  upper  mantle  response  is  nearly  independent 
of  frequency  and  the  approximation  (3.3)  is  not  bad.  The 
instrument  response  was  taken  to  be  the  LRSM  short  period 
nominal  response  for  all  the  computations. 

Seismograms  were  computed  at  depths  ranging  from  0.3 
to  3.0  km.  The  P-pP  lag  times  then  varied  between  0.13  and 
1.31  seconds.  The  suite  of  seismograms  was  computed  for 
three  cases  with  differing  dominant  frequencies.  These  are: 
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TABLE  4 

SOURCE  CRUSTAL  MODEL  FOR  DEPTH  OF  BURIAL  STUDY 


Depth  (km) 

Thickness  (km) 

a (km/sec) 

6 (km/sec) 

p (qm/cm3) 

3.1 

3.1 

4.3 

2.4 

2.6 

7.0 

3.9 

4.7 

2.6 

2.6 

13.0 

6.0 

5.4 

2.7 

2.7 

20.0 

7.0 

6.0 

3.5 

2.8 
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1. 

Yield  = 300 

kt; 

T/Q  = 1.05, 

2. 

Yield  = 300 

* 

r+ 

T/Q  = 0.7, 

3. 

Yield  = 10 

kt? 

T/Q  = 1.05. 

Other  than  the  yield  and  T/Q,  the  parameters  for  the  three 
sets  of  seismograms  are  identical  in  all  respects. 

The  synthetic  seismograms  for  the  three  cases  listed 
above  are  shown  in  Figures  4.1  through  4.3.  The  P-pP  lag 
is  indicated  on  each  record.  To  provide  a standard  of  compari- 
son we  also  computed  the  direct  P wave  seismogram  for  each 
case  and  this  is  shown  at  the  top  of  each  figure.  The  direct 
P wave  computation  was  done  by  taking  the  first  layer  of  the 
crustal  structure  of  Table  4 to  be  a halfspace  rather  than 
having  a free  surface. 

Standard  measurements  of  the  b and  d phases  were  made. 
The  cycle  for  the  d phase  measurement  is  indicated  on  each 
record  by  a bar  at  the  zero  crossing.  Since  these  are 
theoretical  seismograms  given  in  digital  form,  the  peak-to- 
peak  amplitude  and  time  (half  the  apparent  period  of  the 
cycle)  can  be  measured  with  extraordinary  accuracy.  In 
fact,  measurements  are  generally  done  automata o^iiy  using 
parabolic  fits  at  each  peak  rather  than  by  hand.  As  the 
P-pP  lag  time  increases , phasing  becomes  apparent  in  the 
maximum  cycle.  For  some  of  these  seismograms  we  continued 
to  use  the  peak-to-peak  amplitude  but  took  the  half  period 
from  a zero-to-zero  measurement  on  the  uncontaminated 
portion  of  the  wavelet.  The  seismograms  where  this  was 
done  are  indicated  by  an  H symbol  above  the  measured  half 
period.  A consistent  determination  of  this  period  is  quite 
important  because  of  the  instrument  correction. 

The  measured  periods  of  the  b and  d phases,  Tfa  and 
T^,  are  given  on  each  seismogram.  We  see  that  the  period  of 
the  b phase  increases  steadily  until  the  lag  is  about  0.53  - 
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Figure  4.2.  Seismograms  identical  to  those  of  Figure  4.1 
except  T/Q  =0.7.  The  yield  is  300  kt. 
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0.61  seconds,  then  decreases  a bit  to  its  final  level.  The 
d phase  period  behaves  more  erratically  as  the  P and  pP 
phases  pull  apart  and  we  change  the  quantity  being  measured. 
Still  these  measurements  are  probably  the  most  consistent 
and  accurate  possible.  We  also  point  out  that  at  the 
greatest  lag  time  the  P and  pP  phases  are  essentially 
separated  as  can  be  seen  by  comparing  to  the  direct  P wave 
seismogram  at  the  top  of  each  figure. 

Our  objective  is  to  construct  curves  showing  how  pP 
changes  the  teleseismic  amplitudes  from  the  values  for  the 
direct  P wave  taken  alone.  As  we  saw  in  Section  III,  it  is 
the  direct  P wave  that  is  most  indicative  of  yield  if  we 
know  the  properties  of  the  emplacement  material.  Therefore, 
we  normalize  the  b and  d amplitudes  from  the  seismograms  of 
Figures  4.1  - 4.3  to  the  values  for  the  direct  P wave  taken 
separately.  The  results  are  plotted  versus  P-pP  lag  time 
in  Figures  4.4  and  4.5. 

For  the  b phase  corrections,  Figure  4.4,  we  see  that 
the  behavior  is  not  markedly  different  for  the  three  cases 
studied.  At  very  short  lag  times  the  b phase  amplitude  is 
reduced  by  factors  ranging  down  to  0.6  and  below.  We  should 
point  out  that  lag  times  less  than,  say,  0.4  seconds  are 
extremely  unlikely  for  the  large  explosions  of  interest. 

For  the  300  kt  explosion  the  b phase  normalized  amplitude 
peaks  at  around  1.2  so  the  amplitude  is  not  significantly 
enhanced  for  any  lag  time.  The  peaking  for  the  10  kt  event 
is  somewhat  higher  but  the  noisy  character  of  the  records 
makes  the  results  a bit  suspect.  In  summary,  we  conclude 
that  for  likely  values  of  the  P-pP  lag,  the  b phase  is 
closely  representative  of  the  direct  P wave. 

The  d phase  correction  curves.  Figure  4.5,  are  quite 
interesting  because  the  correction  can  be  significant.  First 
we  note  that  except  for  lag  times  in  the  0.6-0. 9 second 
range  the  curves  are  pretty  much  the  same  for  all  three  cases. 
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In  this  lag  time  range  the  P and  pP  phases  separate  and  the 
apparent  period  measurement  cannot  be  made  consistently. 

The  character  of  the  normalized  amplitude  curves,  especially 
the  sharp  discontinuities,  is  probably  more  due  to  rapid 
changes  in  the  instrument  correction  factor  than  to  changes 
in  the  amplitude  of  the  actual  ground  motion.  Still,  we  are 
following  conventional  procedure. 

The  peak  of  the  d amplitude  correction  curve  occurs 
at  about  0.6-0. 7 seconds  for  the  300  kt  events.  At  the  peak 
the  amplitude  is  more  than  1.8  times  as  large  as  that  from  the 
direct  P wave  itself.  More  important,  the  correction  curve 
falls  away  rapidly  on  either  side  of  the  peak  so  events  with 
lags  near  these  values  can  appear  considerably  larger  than 
those  with  lags  a few  tenths  of  a second  different. 

For  the  10  kt  event  the  d amplitude  correction  curve 
peaks  at  an  extraordinarily  high  value,  nearly  2.4.  The 
high  value  is  due  to  the  long  apparent  period  of  the  d phase 
and  commensurate  instrument  correction  factor.  Is  this  an 
aberation?  Perhaps  not.  The  only  different  between  cases 
1 and  3 is  that  the  source  spectrum  is  shifted  to  higher 
frequencies  by  a factor  of  3.1. 

In  the  following  section  the  practical  consequences 
of  the  depth  correction  factors  presented  here  will  be 
explored  by  applying  them  to  actual  data. 

SUMMARY 

If  the  d amplitude  is  used  to  obtain  m^,  then  both 
the  pP  interval  time  and  the  instrument  response  characteris- 
tics must  be  available  before  this  measurement  can  be  used  to 
obtain  the  amplitude  of  the  direct  P wave. 
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V.  MAGNITUDE- YIELD  SCALING  INCLUDING  CORRECTIONS  FOR 

DEPTH  OF  BURIAL 


In  this  section  we  apply  the  theoretical  P-pP  lag 
time  corrections  developed  in  Section  IV  to  a small  set 
of  actual  data.  The  results  are  therefore  quite  preliminary, 
but  intriguing. 

The  data  is  that  compiled  by  Alewine,  et  al . (1976) 

for  a group  of  Pahute  Mesa  events.  This  data  is  presented 
in  the  form  of  magnitude  (m^)  estimates  based  on  the  maxi- 
mum amplitude  (d  phase)  within  the  first  three  cycles  of  the 
P wavetrain  recorded  at  teleseismic  distances.  The  magni- 
tude estimates  were  determined  from  measurements  at  a number 
of  stations,  carefully  averaged  to  remove  sources  of  error 
or  bias.  From  the  data  set  we ‘select  a number  of  large 
yield  events  detonated  below  the  water  table  at  Pahute  Mesa, 
NTS.  Thus,  it  would  seem  that  the  normalized  d amplitude 
curve  for  300  kt,  T/Q  - 1.05  from  Figure  4.5  is  quite  ap- 
propriate and  will  be  used  to  correct  the  observed  m^  values. 

In  order  to  apply  the  corrections,  estimates  of  the 
P-pP  lag  times  are  required  for  the  events  under  considera- 
tion. Several  sources  were  used  to  estimate  this  quantity. 

First,  we  need  delay  times  from  published  reports  by 
Kulh^nek  (1971)  and  Cohen  (1975) . These  authors  employed 
the  spectral  null  technique  in  order  to  determine  P-pP  delay 
times  for  several  Pahute  Mesa  events.  Using  these  delay 
times  and  the  amplitude  correction  curve  of  Figure  4.5,  the 
magnitude  estimates  from  Alewine,  et  al . were  modified. 

The  corrected  m^  values  are  plotted  versus  explosion  yield 
in  Figure  5.1  as  closed  circles.  The  dashed  open  circles 
are  the  original  (uncorrected)  estimates  from  Alewine,  et  al . 

Least  squares  fits  were  made  to  each  set  of  data  points 
in  Figure  5.1.  The  results  are  as  follows: 
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Original  Estimates 

Adjusted  for  Reported 
Spectral  Null  pP-P 
Delay  Times 


Yield  (kt) 


Original  and  revised  (corrected)  magnitude  esti- 
mates versus  yield  for  a population  of  Pahute 
Mesa  explosions.  The  P-pP  lag  times  for  the 
corrections  are  based  on  spectral  null  and  de- 
convolution  techniques  and  are  given  by 
Kulhanek  (1971)  , Cohen  (1975)  and  Frasier  (1972) 


Original  Estimates 

@ Adjusted  for  pP-P  Delay 
Tildes  Eased  on  CEP 
Structures 
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Yield  (kt) 


Original  and  revised  (corrected)  magnitude  esti- 
mates versus  yield  for  a population  of  Pahute 
Mesa  explosions.  The  P-pP  lag  times  for  the 
corrections  are  based  on  the  sonic  velocities  of 
the  overburden  material. 
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1.  Adjusted  ir^'s,  = 3.45  + 0.98  log  Y, 

2.  Original  m^'s,  = 4.35  + 0.66  log  Y. 

The  most  significant  aspect  of  the  data  in  Figure  5.1  is 

the  fact  that  the  slope  of  the  straight  line  fit  to  the  ad- 

justed m^  values  is  quite  close  to  unity  (0.98),  as  com- 
pared to  the  much  smaller  slope  for  the  original  m^  estimates. 
Values  of  m^  - log  Y slopes  between  0.75  and  0.95  are  quite 
common  in  the  literature  for  what  we  have  here  termed  un- 
corrected m^  estimates.  It  seems  that  the  pP  effect  is  a 
strong  contributing  factor  to  this  deviation  from  unit 
slope . 

A second  set  of  pP-P  delay  times  was  estimated  for  a 
larger  (13  events)  population  of  Pahute  Mesa  explosions. 

These  delay  times  were  based  upon  best  estimates  of  the 
vertical  two-way  travel  time  between  the  explosion  working 
point  and  the  free  surface.  Detailed  crustal  structures 
from  well  log  data  were  available  for  two  of  the  explosions 
in  this  population.  The  remaining  11  explosions  were  deto- 
nated in  the  vicinity  of  one  or  the  other  of  these  two  shots 
and  the  corresponding  structures  used  to  estimate  travel 
times.  In  general,  the  delay  times  based  on  these  crustal 
structures  proved  to  be  significantly  shorter  (^  0.2  seconds) 
than  those  based  on  the  spectral  null  techniques  of 
Kulhanek  (1971)  and  Cohen  (1975)  or  the  deconvolution  pro- 
cedure described  by  Frasier  (1972)  . Using  these  estimates 
or  delay  times  and  the  amplitude  correction  curve  in  Figure 
4.5,  magnitude  estimates  given  by  Alewine,  et  al . (1976) 

were  adjusted  and  plotted  versus  yield  in  Figure  5.2. 

The  effect  of  the  P-pP  lag  time  correction  on  this 
population  of  events  is  significantly  different  from  that 
in  Figure  5.1.  In  this  case,  the  corrections  are  all  in  the 
same  direction  and  least  square  fits  to  original  and  corrected 
data  give  about  the  same  results: 
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Adjusted  m^'s,  = 3.45  + 0.90  log  Y, 

2.  Original  ra^'s,  m^  = 3.61  + 0.91  log  Y. 

What  is  the  meaning  of  the  contradictory  results  of 
Figures  5.1  and  5.2?  First  consider  the  results  of  Figure 
5.1  which  were  based  on  estimates  of  the  P-pP  lag  time  taken 
directly  from  recorded  seismograms.  While  these  may  not  be 
terribly  accurate  in  any  particular  example,  a number  of 
authors  have  obtained  about  the  same  results.  Further,  these 
lag  times  are  not  far  different  from  those  inferred  from 
direct  measurements  of  the  source-to-surface  travel  time 
when  the  explosion  is  detonated  (Springer,  1974). 

lefer  to  the  normalized  d amplitude  correction  curve 
of  Figure  4.5.  The  reason  for  the  change  in  slope  that 
occurs  in  the  data  of  Figure  5.1  when  the  P-pP  lag  correc- 
tion is  applied  is  quite  simple.  The  larger  yield  explosions 
fall  on  the  right  side  of  the  peak  where  the  correction  is 
small  while  the  smaller  yield  explosions  fall  much  nearer 
the  peak  on  the  left  side  where  the  correction  is  relatively 
larger. 

Estimates  of  the  P-pP  lag  based  on  the  sonic  velocity 
properties  of  the  overburden  rock  are  uniformly  smaller  than 
the  estimates  based  on  actual  observations  of  the  explosions. 
It's  not  entirely  clear  why  this  occurs,  but  it  seems  to  be 
a fact.  If  we  use  the  much  lower  sonic  estimates  of  the  P- 
pP  lag,  all  the  corrections  to  the  data  of  Figure  5.2  fall 
on  the  left  side  of  the  peak  in  Figure  4.5.  All  the  events 
are  then  corrected  downward  by  a fairly  uniform  increment 
and  the  m^-logY  slope  is  not  changed. 

As  we  said  at  the  beginning  of  the  section,  these 
results  are  preliminary,  but  quite  intriguing.  The  corrections 
have  been  applied  to  only  a few  events  in  a single  test  re- 
gion. The  study  should  be  broadened  to  include  many  more 
events  in  different  test  regions.  Perhaps  this  will  make  a 
significant  contribution  to  resolving  some  of  the  anomalies 
in  the  data. 
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SUMMARY 


What  is  the  best  way  to  obtain  the  pP  interval  time? 

If  we  use  spectral  nulls  and/or  deconvolve  the  recorded 
ground  motion  then  we  no  longer  have  a predictive  capability. 
Furthermore,  this  technique  is  subject  to  evasion  possibili- 
ties. 

We  recommend  that  our  shot  data  be  used  to  determine 
the  in  situ  P wave  velocity  measurements  which  best  agree 
with  the  analysis  of  recorded  teleseismic  ground  motion. 

If  we  find  that  we  need  the  source-to-surface  travel  time 
for  the  calibration  shots  we  should  ask  for  it. 
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VI.  TWO  "ANOMALOUS"  EVENTS;  FAULTLESS  AND  A FRENCH 
EXPLOSION  AT  THE  HOGGAR  TEST  SITE 


FAULTLESS  is  considered  to  be  anomalous  because  it 
falls  considerably  above  the  empirical  m^-log  Y curve  based 
on  events  in  similar  source  materials.  Alewine,  et  al. 

(1976)  give  m^  = 6.51  for  FAULTLESS.  Springer  (1974)  esti- 
mated the  P-pP  lag  time  for  this  event  to  be  0.8  seconds 
using  the  measured  arrival  time  of  the  explosion  shock  wave 
at  surface  zero.  Frasier  (1972)  gives  0.90  seconds  and 
Kulhanek  (1971)  gives  1.02  seconds  for  the  same  quantity  us- 
ing teleseismic  data.  If  we  use  Springer's  value  of  0.80  and 
the  correction  curve  for  Y = 300,  T/Q  = 1.05  from  Figure 
4.5,  the  revised  magnitude  estimate  is  6.31,  lower  than  the 
original  estimate  by  0.2  units. 

If  the  revised  magnitude  for  FAULTLESS  is  plotted  in 
Figure  5.1,  the  anomalous  magnitude-yield  scaling  for  this 
event  practically  disappears.  The  revised  magnitude  plots 
within  0.1  magnitude  unit  of  the  least  square  fit  (either 
m^  = 3.09  + 1.1  log  Y or  m^  = 3.45  + 0.98  log  Y)  to  the 
corrected  magnitude  data  shown  in  Figure  5.1.  Therefore, 
if  we  apply  the  P-pP  lag  correction  in  a consistent  way  to 
Faultless  and  a suite  of  Pahute  Mesa  events  of  similar  yield 
and  in  similar  emplacement  materials,  the  FAULTLESS  event 
no  longer  appears  to  be  anomalous. 

The  French  explosions  at  the  Hoggar  site  in  the 
Sahara  are  considered  to  be  anomalous  because  they  are  es- 
sentially undistinguishable  from  United  States  granite  explo- 
sions (e.g. , PILEDRIVER)  or  an  m^-log  Y plot.  However, 
the  French  report  that  cavity  radii  in  the  Hoggar  granite 
are  40  percent  smaller  then  U.S.  granite  cavity  radii,  i.e., 
PILEDRIVER  and  HARDHAT . The  source  functions  for  events  in 
NTS  and  Hoggar  granite  were  presented  in  Section  II.  We 
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saw  that  the  Hoggar  source  function  is  much  smaller;  the 

A 

ct  | H7  ( f ) | is  down  by  a factor  of  about  5 or  6 in  the  fre- 
quency  range  of  interest.  Further,  we  know  that  the  NTS 
granite  source  gives  a good  match  to  both  near  field  and 
teleseismic  data  for  the  PILEDRIVER  event  (Bache,  et  al . , 
1975) . How  then  can  we  explain  the  fact  that  the  m^  is  about 
the  same  for  explosions  of  the  same  yield  in  the  two  grani- 
tic materials? 

The  first  step  to  be  taken  is  to  see  how  much  we 
really  do  underestimate  the  Hoggar  body  waves  when  all  fac- 
tors are  taken  into  consideration.  We  have  available  a re- 
cording of  a February  27,  1965,  Hoggar  explosion  from  the 
LRSM  station  HNME.  We  also  have  an  estimate  of  the  yield 
of  this  explosion.  We  then  proceed  to  synthesize  a seismo- 
gram at  HNME  using  the  Hoggar  granite  source  (202)  , shown 
in  Figure  2.1. 

The  factors  we  do  not  know  are  the  crustal  structure 
at  the  source  (important  only  for  the  P-pP  lag  time) , the 
upper  mantle  response  and  the  best  value  of  T/Q.  We  do  not 
really  know  the  crustal  structure  at  the  receiver  or  the 
instrument  response  either  but  we  choose  to  use  the  same 
models  as  in  our  other  work.  The  instrument  response  is 
taken  to  be  the  LRSM  nominal.  For  the  crust  at  the  receiver 
we  used  the  same  transparent  crust  used  in  the  synthetic 
seismogram  calculations  of  Section  IV.  This  is  tabluated  in 
Table  5.  The  elementary  source  crustal  model  used  also 
appears  in  this  table. 

The  upper  mantle  response  is  required  at  a distance 
(z  62°)  outside  our  previous  experience.  We  also  wish  to 
synthesize  records  at  another  LRSM  station,  RKON,  which  is 
even  farther  (=  77°).  One  of  the  best  available  upper 
mantle  models  is  C2,  recently  published  by  Anderson  and 
Hart  (1976) . The  upper  mantle  response  at  the  HNME  and 
RKON  distances  was  computed  using  our  generalized  ray 
theory  programs  and  this  model.  We  show  in  the  next  section 


39 


I 


TABLE  5 

CRUSTAL  STRUCTURES  FOR  SYNTHETIC  SEISMOGRAMS  FOR 
EXPLOSIONS  AT  THE  HOGGAR  TEST  SITE 


Source  Structure 

Depth  (km) 

a (km/sec) 

B (km/sec) 

p (gm/cm3) 

2.40 

5.33 

3.78 

2.67 

5.20 

5.80 

3.45 

2.67 

21.20 

6.00 

3.50 

2.80 

Receiver  Structure 

2.58 

3.67 

2.31 

2.40 

4.84 

5.42 

3.27 

2.60 

11.61 

5.80 

3.45  . 

2.60 

20.00 

6.00 

3.50 

2.80 
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that  in  this  range  the  variation  in  amplitude  from  model 
to  model  is  likely  to  be  less  than  a factor  of  two. 

Having  disposed  of  all  the  "elastic"  parameters, 
there  are  two  "inelastic"  parameters  to  be  specified,  T/Q 
and  the  P-pP  lag  time.  The  latter  is  classified  as  inelas- 
tic because  we  believe  that  inelastic  effects  are  responsible 
for  the  lag  time  being  considerably  longer  than  that  obtained 
from  sonic  velocities  of  the  actual  overburden  material. 

These  two  quantities  shape  the  seismogram  and,  to  a consider- 
able degree,  determine  its  amplitude. 

First,  consider  T/Q.  It  is  widely  believed  that  NTS 
is  underlain  by  a region  of  anomalously  low  Q.  Therefore, 
we  would  expect  to  require  a Hoggar  value  of  T/Q  lower  than 
the  0.95  - 1.05  values  we  have  found  appropriate  for  NTS 
events.  The  frequency  content  of  the  observed  seismogram 
also  strongly  constrains  T/Q. 

For  the  P-pP  lag  time  the  only  constraint  is  the 
shape  of  the  recorded  seismogram  and,  perhaps,  past  experience 
with  events  in  similar  materials.  We  synthesized  seismograms 
for  a number  of  lag  time  values  and  finally  settled  0.80 
seconds  as  giving  the  best  match  to  the  seismogram.  This 
lag  time  is  consistent  with  estimates  for  the  lag  time  for 
PILEDRIVER  made  from  teleseismic  data.  Since  we  don't 
know  the  burial  depth  of  the  Hoggar  event,  we  can  only 
conjecture  regarding  its  depth  relative  to  Piledriver. 

Synthetic  seismograms  were  computed  at  HNME  for  two 
values  of  T/Q;  0.7  and  0.5.  The  two  seismograms  are  shown 
in  Figure  6.1  together  with  the  observed  record.  The  T/Q  = 

0.7  seismogram  matches  the  frequency  content  of  the  first 
few  seconds  of  the  observed  record  quite  well.  However, 
the  amplitude  is  clearly  too  small.  The  T/Q  = 0.5  record 
is  useful  because  it  illustrates  the  fact  that  this 
quantity  is  narrowly  constrained  by  the  frequency  content 
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Comparison  of  synthetic  (bottom  two)  and  ob- 
served (top)  short  period  vertical  seismograms 
at  HNME  for  the  27  February,  1965,  Hoggar 
explosion.  The  time  and  amplitude  scales  are 
the  same  on  all  seismograms.  The  amplitude 
in  microns  at  1 Hz  is  indicated  at  the  left 
on  the  synthetic  records. 


of  the  observations.  Even  though  the  amplitude  of  the 
T/Q  = 0.5  record  more  closely  matches  the  observed  ampli- 
tude, the  dominant  frequencies  are  clearly  much  too  high. 

This  fact  cannot  be  changed  by  adjusting  the  P-pP  lag  time. 

The  synthetic  seismogram  with  a P-pP  lag  time  of  0.8 
seconds  and  T/Q  = 0.7  gives  an  excellent  match  to  the  fre- 
quency content  of  the  HNME  observation.  What  can  we  say 
about  the  amplitude  and,  in  particular,  m^?  There  are  two 
phases  from  which  m^  can  be  measured  on  the  observation. 

There  are  the  b phase  and  the  d phase  which  is  indicated  by 
a bar  on  the  records  of  Figure  6.1.  On  the  synthetic  seismo- 
gram both  the  b and  d phases  give  the  same  m^,  which  is  5.57. 

For  the  observed  HNME  seismogram  we  compute  m^  - 5.78 
for  the  b phase  and  m^  = 5.93  for  the  d phase.  Based  on 
empirical  m^-log  Y plots  on  which  this  event  appears,  the 
average  m^  seems  to  be  about  5.70.  We  will  therefore  assume 
that  most  of  the  m^  peaks  were  from  the  b phase.  With  this 
assumption  the  HNME  record  is  typical  with  m^  being  a little 
less  than  0.1  units  above  the  average. 

In  summary,  the  theoretical  m^  for  HNME  is  about  0.2 
m^  units  below  the  observed  m^  for  this  station.  The  dis- 
crepency  is  far  less  than  the  factor  of  five  or  six  expected 
from  comparison  of  the  source  functions.  Further,  the  synthe- 
tic m^  is  only  0.13  units  lower  than  the  average  m^  for  the 
Hoggar  event.  In  short,  our  theoretical  amplitudes  match 
the  observations  to  a factor  of  considerably  less  than  two. 

Seismograms  appropriate  for  the  RKON  station  were  also 
computed  with  T/Q  = 0.7  and  0.5  and  these  are  shown  in  Figure 
6.2.  For  the  T/Q  = 0.7  case,  the  theoretical  m^  is  5.43  for 
RKON.  This  is  somewhat  lower  than  that  of  HNME  but  is  still 
within  0.3  m^  units  (or  a factor  of  two  in  amplitude)  of  the 
average  for  the  event. 
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We  synthesize  the  seismograms  in  Figures  6.1  and  6.2 
to  see  how  much  we  actually  underestimate  the  Hoggar  explo- 
sion body  wave  amplitudes.  The  answer  is,  not  very  much. 

As  with  FAULTLESS,  we  see  that  the  Hoggar  events  are  probably 
not  anomalous  to  any  serious  degree. 

SUMMARY 

The  FAULTLESS  anomaly  can  be  explained  by  a pP  correc- 
tion. The  Hoggar  "anomaly"  can  be  explained  by  using  a T/Q 
of  0.7. 


VII.  EFFECT  OF  THE  UPPER  MANTLE  ELASTIC  PROPERTIES 
ON  SHORT  PERIOD  SEISMOGRAMS 


In  this  section  we  address  the  question,  how  much 
can  the  amplitude  and  waveform  of  teleseisms  be  influenced 
by  reasonable  changes  in  the  upper  mantle  velocity  profile? 

We  restrict  our  attention  to  ranges  between  20  and  36°  and 
show  the  effect  by  computing  synthetic  seismograms  with  all 
factors  held  constant  except  the  upper  mantle  velocity  pro- 
file. The  example  chosen  is  the  MAST  explosion  in  rhyolite 
which  was  studied  in  some  detail  by  Barker,  et  al . (1976). 

The  upper  mantle  models  studied  are  perturbations  of 
models  HWNE  (Helmberger  and  Wiggins,  1971)  and  HWA  (Wiggins 
and  Helmberger,  1973) . The  perturbations  were  not  allowed 
outside  the  envelope  of  all  possible  models  consistent  with 
travel  time,  ray  parameter  and  amplitude  data  for  the  western 
United  States,  as  determined  by  Wiggins,  et  al . (1973).  The 
changes  from  the  starting  models  (which  are  fairly  different 
from  each  other)  were  at  most  0.3  km/sec. 

The  synthetic  MAST  seismograms  were  all  computed  with 
the  MAST  source  (Figure  2.1),  T/Q  = 1.05  and  the  instrument 
response  was  the  LRSM  nominal.  Seismograms  were  computed 
for  many  models  but  only  four  examples  which  represent  the 
extremes  in  their  effect  are  discussed  here.  These  models 
are  denoted  HWA-2,  HWA-3B,  HWA-6D  and  HWNE-3  and  their 
velocity-depth  profiles  are  shown  in  Figure  7.1. 

Synthetic  seismograms  appropriate  for  the  epicentral 
distance  from  NTS  to  the  SDCS  stations  RKON , WH2YK  and  HNME 
are  shown  in  Figures  7.2  through  7.4.  These  seismograms 
represent  the  extremes  of  the  variation  observed  at  these 
stations  when  perturbing  the  models  under  the  constraints 
mentioned  above.  The  signals  at  RKON  and  WH2YK  are  affected 
primarily  by  the  structure  100  km  above  and  below  the  400  km 
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discontinuity.  The  HNME  seismograms  are  primarily  sensitive 
to  the  region  below  the  600  km  discontinuity. 

The  b and  d amplitudes  were  measured  on  each  seismo- 
gram and  are  listed  on  Figures  7.2  through  7.4.  At  the 
stations  within  the  triplication  range  (RKON  and  WH2YK)  the 
ratio  of  the  maximum  to  minimum  amplitudes  is  about  a fac- 
tor of  2-2.5,  though  the  b phase  for  HWA-3B  at  RKON  is 
exceptionally  small.  Beyond  the  triplications  at  HNME,  the 
variation  is  less,  a factor  of  1.5-1. 6. 

The  small  variation  in  amplitude  at  the  station  beyond 
the  triplications  is  reassuring.  It  is,  to  be  sure,  based 
on  a limited  number  of  model  perturbations,  but  many  of  these 
perturbations  were  designed  to  change  the  amplitude  at  the 
HNME  range.  Any  model  used  must  be  consistent  with  the 
available  travel  time  data,  etc.,  and  these  data  severely 
constrain  the  velocity  profile.  If  we  assume  that  parallel 
layered  models  adequately  represent  the  upper  mantle  for  our 
purposes,  we  can  be  reasonably  certain  that  errors  in  our 
knowledge  of  the  velocity-depth  profile  lead  to  errors  in 
our  computed  amplitudes  that  are  on  the  order  of  a factor 
of  1.5  to  1.6. 

If  we  are  within  the  triplication  range,  what  are  the 
possible  variations  that  occur  with  small  changes  in  range? 
Once  again,  we  provide  at  least  a partial  answer  to  the 
question  by  showing  a typical  example.  Seismograms  for  the 
MAST  event  were  again  computed  for  a series  of  epicentral 
distances  near  that  for  RKON.  The  model  is  HWNE-3.  The 
seismograms  together  with  the  b and  d amplitudes  are  shown 
in  Figure  7.5.  Over  the  67  km  range  represented  the  rela- 
tive amplitudes  and  arrival  times  of  the  incoming  rays 
vary  considerably  as  evidenced  by  the  changing  shape  of  the 
seismograms.  The  ratio  of  the  maximum  to  minimum  b and  d 
amplitudes  is  1.5  and  1.8,  respectively. 
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